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ABSTRACT 

The degree to which interstellar grains align with respect to the interstellar mag- 
netic field depends on disaligning as well as aligning mechanisms. For decades, it was 
assumed that disalignment was due primarily to the random angular impulses a grain 
receives when colliding with gas-phase atoms. Recently, a new disalignment mechanism 
has been considered, which may be very potent for a grain that has a time-varying 
electric dipole moment and drifts across the magnetic field. We provide quantitative 
estimates of the disalignment times for silicate grains with size > 0.1/im. These appear 
to be shorter than the time-scale for alignment by radiative torques, unless the grains 
contain superparamagnetic inclusions. 
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1 INTRODUCTION 

Observations of starlight polarization have revealed that 
some interstellar dust grains are non-spherical and aligned. 
The degree of alignment, and hence the polarization, de- 
pends on both aligning processes (e.g., radiative torques and 
paramagnetic dissipation) and disaligning processes (e.g., 
random torques arising from collisions with gas atoms). See 
Whittet (2004) for a review of polarization observations and 
Lazarian (2003), Roberge (2004), and Lazarian (2007) for 
reviews of alignment theory. 

Recently, Weingartner (2006, hereafter W06) proposed 
an alternative disalignment mechanism for a grain that has a 
time-varying electric dipole moment p and drifts across the 
interstellar magnetic field. The potency of this mechanism is 
sensitive to the magnitude and time-scale of fluctuations in 
p. W06 considered highly simplified models for the fluctu- 
ating electric dipole moment. Here, we examine this process 
and the implications for disalignment in greater detail. We 
will consider relatively large (size >0.1 fim) silicate grains, 
since the 9.7 and 20pm features exhibit polarization (e.g., 
Martin & Whittet 1990; Smith et al. 2000; Wright et al. 
2002) and the wavelength dependence of the observed po- 
larization implies that relatively small grains are not efficient 
polarizers (Kim & Martin 1995). 

In Sj2[ we review the main elements of disalignment asso- 
ciated with fluctuations in p. Next, we introduce models for 
the transport of charge to and within silicate grains (33]). 
We describe simulations of the fluctuating dipole moment 
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and associated disalignment in 
and conclusions in 36] 



]and present results in 



2 DISALIGNMENT ASSOCIATED WITH 
TIME- VARYING ELECTRIC DIPOLE 
MOMENTS 

When a gas atom collides with a grain, it imparts an angular 
impulse to the grain. If no other mechanisms excite rotation, 
then the energy in rotation about any axis is ~ |/cBT gas , 
where fee is Boltzmann's constant and T gas is the gas tem- 
perature. Such motion is called 'thermal rotation'. The ther- 
mal rotation rate for a sphere with radius a is given by 
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where p is the density of the grain material. In general, grains 
are subjected to additional torques that may drive them to 
suprathermal rotation, with angular speed uj > ujt (Purcell 
1975, 1979; Draine & Lazarian 1998). For thermally rotating 
grains, the random collisional impulses constitute an impor- 
tant disalignment mechanism. 

A spinning grain with non-zero electric charge acquires 
a magnetic dipole moment \i || u (Martin 1971). Dolginov 
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& Mytrophanov (1976) showed that the Barnett effect (i.e., 
the tendency for a spinning paramagnetic solid to acquire 
a magnetization parallel or anti-parallel to ui) can provide 
a much larger moment. Specifically, the Barnett magnetic 
moment //Bar = XowV/7 s , where xo is the static magnetic 
susceptibility, 7 9 is the gyromagnetic ratio of the microscopic 
magnetic dipoles that are responsible for the grain's para- 
magnetism, and V is the grain volume. The susceptibility 
depends on the number density of paramagnetic ions or nu- 
clei in the grain material and is thus rather uncertain. We 
will adopt xo ~ 5 x 10~ 3 (T d /15 K)~\ where T d is the grain 
temperature (Draine 1996; W06). With this estimate, the 
Barnett magnetic moment for silicate grains is 



|MBar|(sil) 




Suppose fj, is parallel or anti-parallel to the grain's an- 
gular momentum vector J: \i = \ijJ. The magnetic torque 
T M = (a X B causes J to precess about the interstellar mag- 
netic field B at rate 




Since p oc ui, Qo is independent of to. The combination of 
an aligning torque (e.g., the radiative torque) and the mag- 
netic torque drives the grain towards rapid precession with 
a constant precession angle # a ii g n (i.e., (9 a ii g n is the angle 
between B and J). A large ensemble of grains will be char- 
acterized by a uniform distribution in precession phase. As 
a result, the observed starlight polarization is either parallel 
or perpendicular to B. (If # a ii gn = and J || &i, the grain 
principal axis of greatest moment of inertia, then the polar- 
ization || B.) Note that alignment of the grain body with 
respect to J is also a necessary condition for polarization. 

If a grain has an electric dipole moment p and drifts 
with velocity v across B, then it experiences a torque T p = 
p X (v X B)/c in addition to the magnetic torque (c is the 
speed of light). If p = pjJ, then the grain precesses about 
an axis tilted at angle 5 = tan -1 |T| relative to B and the 
precession rate is increased by the factor (1 + T 2 ) 1 / 2 , where 

y _ Pjv±_ (4) 

with v± the component of v transverse to B (W06). 

If (j, and p are not parallel or anti-parallel to J, then 
the magnetic and electric torques must be averaged over 
the extremely rapid grain rotation. The resulting dynamics 
is identical to that for which J , fi, and p all lie along &i, 
except with the following substitutions in equation Q: 

<«-(£)£ ■ < 5 > 

PJ = ±(p- a,i) fi(q) , (6) 

with Ii the moment of inertia along at (W06). Both \ij 
and pj depend on the grain's rotational state through the 



parameter q = 2IiE/J 2 (E is the rotational energy). In 
equation (|6}, the + (-) sign is selected when J -en > (< 0) 
and the factor fi(q) is given in eq. 9 of W06. (The choice 
of i is also discussed following eq. 9 in W06.) If the Barnett 
effect is responsible for the magnetic dipole moment, then 

O/M = XoV/lg- 

If T is constant in time, then the only consequence of 
the electric dipole is to tilt the precession axis relative to the 
magnetic field direction. However, T can vary on relatively 
short time-scales. W06 discussed two sources of variation: 1. 
Upon each discrete charging event (e.g., the capture of an 
electron from the gas or photoejection of an electron), pj 
changes. 2. When the grain's rotational state (q and/or flip 
state) varies, p.j / p.j varies (eqs. [5] and [6]). 

Two processes can yield rapid variations in the grain 
rotational state: 1. Thermal fluctuations, in which energy 
is exchanged between grain rotation and vibrational modes 
(Lazarian 1994; Lazarian & Roberge 1997; Lazarian & 
Draine 1997, 1999a, 1999b; Weingartner 2009). 2. Collisions 
with gas-phase atoms, which can stick to, reflect from, or 
evaporate from the surface, perhaps after forming a molecule 
(Hoang & Lazarian 2009). The efficacy of both of these 
mechanisms drops off dramatically as the grain rotation be- 
comes suprathermal. 

In this paper, we will only consider variations in pj as- 
sociated with discrete charging events. We also will assume 
q = 1 in equation <(Sj and fi(q) = 1 in equation ([6]), which 
are good approximations for suprathermally rotating grains. 
Thus, the analysis presented here is not complete for ther- 
mally rotating grains. 

Since the charging processes are stochastic processes, T 
varies stochastically, yielding random variations in the pre- 
cession axis. Each time the precession axis changes direc- 
tion, the precession angle changes. When these events occur 
at random precession phases, 6* a ii gn varies stochastically. In 
other words, the grain experiences disalignment. 

W06 considered a simple scenario in which T has con- 
stant magnitude but stochastically reverses sign, on time- 
scale Tflip, finding the following approximations for the dis- 
alignment time-scale when rmp is short or long compared 
with the precession time-scale: 

r dis ~ t- 2 in r 2 r M i , if r flip « |fi r l (i+r 2 r i/2 (7) 

r dis ~ (1 + T- 2 )T flip , if rflip > Ifior^l + T 2 )" 1 / 2 . (8) 

In the following section, we will consider more detailed mod- 
els for the fluctuating electric dipole moment. 



3 GRAIN CHARGING MODELS 

A grain charging model that follows the evolution of the 
electric dipole moment p must treat both the processes that 
deliver charge to the grain and those that transport charge 
within the grain. In the cold, neutral, interstellar medium, 
the dominant charge delivery mechanisms are starlight- 
induced photoelectric emission and sticking collisions of gas- 
phase electrons. 
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3.1 Idealizations for Charge Transport Within a 
Grain 

Bulk, neutral silicates are good insulators, with a full va- 
lence band and empty conduction band. Observations of 
the 9.7 /im band profile indicate that interstellar silicates 
are predominantly amorphous (Li & Draine 2001; Kemper, 
Vriend, k, Tielens 2004; Li, Zhao, & Li 2007). In amorphous 
materials, localized energy states ('traps') appear in the tails 
of the conduction and valence bands. For any realistic inter- 
stellar grain, there are also localized states associated with 
impurity atoms. Electrons and holes can hop from site to 
site with assistance from a phonon (e.g., Mott & Davis 1971; 
Blaise 2001), so no grain is perfectly insulating. 

The rate at which an electron hops from site i to site j 
is typically approximated as 

-Rhop = ^ P h exp(-2nj/d^) exp(-Wi^j/k B T d ) (9) 

(Ambegaokar et al. 1971; Mady et al. 2007), where v p h ~ 
10 13 s _1 is the phonon frequency (Brucato et al. 2002), ry 
is the distance between sites i and j, d^ is the electron local- 
ization length, Wi-tj = max[i?j — Ei,0], Ei is the electron 
energy when localized at site i, and Td is the dust tempera- 
ture. 

A completely rigorous treatment of the grain electric 
dipole moment would include following the charges as they 
hop among traps. However, this approach is not feasible. 
First, the quantities appearing in equation ((9j), namely d^> 
and the trap energy distribution, are poorly known. Second, 
even for tight binding at traps (e.g., d^ « 2 A), there are 
typically numerous neighboring traps for which the hopping 
time is orders of magnitude smaller than the time between 
discrete charging events (which itself is orders of magnitude 
smaller than the disalignment time). This is true even when 
a charge is well-localized within the vicinity of a particu- 
larly deep trap. Given the large disparity in time-scales, the 
CPU time for a simulation that follows hopping in detail is 
prohibitive. 

Fortunately, a few simple, plausible idealizations are 
available and do not strain computational resources. We will 
consider the following 4 models: 

1. A perfect insulator. Each time a charge arrives at the 
grain (either an electron from the gas or a hole left following 
photoemission), it remains at its arrival point forever. The 
full hopping model simplifies to this case when — * 0, if the 
typical distance between traps is much less than the grain 
size. This idealization is also reasonable if (a) the typical 
distance between deep traps is much less than the grain size 
and (b) the deep traps effectively retain charges in their 
immediate vicinity. That is, a charge is unlikely to leave the 
'sphere of influence' of a deep trap before recombining. 

2. A perfect conductor. The excess charge on the grain is 
completely delocalized. For a homogeneous, spherical grain, 
the electric dipole moment p vanishes in this case. For non- 
spherical shapes, p oc Z, the net charge on the grain (in units 
of the proton charge). This model is probably not suitable 
for interstellar grains, since we expect Z to be less than the 
total number of deep traps in the grain. Still, it is useful 
to consider this case, to constrain the range of possible out- 
comes. 

3. A conducting grain with deep traps. Some number 
of deep traps are located at random positions within the 



grain. When a charge arrives, it immediately moves to the 
nearest available trap (either occupying it or recombining 
with a resident charge of the opposite sign). This model 
approaches case (1) as the number density of deep traps 



increases. 



4. A partially conducting grain with deep traps. Same 
as (3), except that a charge executes a random walk through 
the grain, with some typical step size and frequency, until it 
comes close to an available deep trap, where it gets stuck. We 
assume that any adsorbates present on the grain surface are 
sufficiently dilute that there is no associated enhancement 
in conductivity along the surface. 



3.2 Collisional Charging 

The trajectories of charged particles in the vicinity of a grain 
with non-vanishing electric dipole moment p differ from 
those for the p — case. The distribution of arrival sites 
on the grain surface is such as to reduce p — \p\. Except for 
model (2) in H3.ll this effect is critical for limiting p. How- 
ever, it is extremely difficult to treat for non-spherical grain 
shapes. Thus, we will always treat the grain as a sphere 
when computing collisional charging rates and the arrival 
sites of colliding particles. For further simplification in these 
calculations, we also neglect the motion of the grain with re- 
spect to the gas. Even though the grain's speed is assumed 
to be roughly the sound speed of the gas, the speed of the 
light electrons is greater by a factor ~ (m p /m e ) 1//2 (m p and 
m e are the proton and electron mass, respectively). Thus, 
we do not expect this assumption to introduce serious er- 
ror for electron collisional charging. In addition, we neglect 
ion collisional charging, which is dominated by photoelectric 
emission. These simplifications are justified in Appendix fAl 
For a grain at rest with respect to the gas, the collisional 
charging rate is given by 



r> 2 /8fcB^as\ 1/2 R 

R — Tva ns — ] R 

\ 7rm 



(10) 



where n is the number density of the colliding particles, 
s is the sticking coefficient (i.e., the probability that the 
particle sticks to the grain following a collision), m is the 
mass of colliding particle, and R accounts for deviations 
of the collision cross section from the geometric cross sec- 
tion. For the relatively large grains under consideration here, 
we adopt s m 1/2 (Weingartner & Draine 2001, hereafter 
WD01). Draine & Sutin (1987) provided expressions for R 
for a charged, conducting sphere, including the polarization 
of the grain by the charged gas-phase particle. The effect of 
polarization decreases with grain size (as long as T gas does 
not approach zero), and can be reasonably neglected when 
a > 0.1 /»m. 

Consider a spherical grain with radius a centered at the 
origin. Approximate the charge distribution within the grain 
as a point charge Q and point dipole pz (p > 0) located at 
the origin. In spherical coordinates (r, 8, <j}), the electric force 
on a point charge q is 



F = ^T* + % (2cos 6»r + sin 0) 
and the potential is 
L t = Qq j IV cos 



(11) 



(12) 
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The equations of motion are 



mr9 2 + mr sin 2 



+ Qq + 2gy cos e 



mr9 - 



-2mr9 + mr sin cos 9 



■ 9 qp sin 9 
t> + ^—^ — 



= — 2mfsin£ 



2mr cos 9 1 



(13) 

(14) 
(15) 



where dots denote differentiation with respect to time and 
m is the mass of the point charge q. 

Employing Hamilton- Jacobi theory, we find the follow- 
ing conserved quantities: 
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(16) 
(17) 
(18) 



as can be verified by direct time differentiation, substituting 
for the second derivatives from equations (|13|l - l|15p . 

Our goals are to find (1) the rate at which incoming 
charged particles strike the grain surface and (2) the distri- 
bution of their arrival angles 9, given Q, p and T g as. First, 
we describe the trajectory of the incoming particle when it 
is still far from the grain (see Fig.[T]). Suppose its velocity is 



= — v(cos 9q z + sin 9q x) 



(19) 



The trajectory is offset from the line x = tan 9o z, which 
passes through the grain center, by impact parameter b; an- 
gle a specifies the displacement of the trajectory from the 
x — z plane. Consider a plane front of incoming particles. 
When the particle whose trajectory passes through the ori- 
gin is located at distance ro from the origin, the coordinates 
of the other particles are 



(r,M)« ( 



r ,9 +cosa — , 



(20) 

(21) 
(22) 
(23) 

Expressing the time derivatives in equations (I16|l - (|18p 
using equations (|21|) (I23|l yields 



ro sin Uo ro , 
and the components of their velocities are 

r m —v 

r9 w v cos a (b/ro) 
rsin9<j> ~ vsina (b/ro). 



p<t> — mvb sin 9o sin a 
P = m 2 v 2 b 2 + 2mqpcos( 
E = mv 2 j2. 



(24) 
(25) 
(26) 



Substituting these results in equations (|17p and (|18p . we find 



V r z J mr \ 



p cos 9o 



(27) 



sin 2 9o sin 2 a 



sin 2 9 



+ — (cos O - cos 0) .(28) 
m 



We are interested in the solution with r < 0, since the par- 
ticle approaches the grain. The choice of the initial sign S 
is more complicated for 9. If 9o = (n), then S = +1 
(5* = —1). Otherwise, equation (|22p yields 9 — vbcosa/r 2 




Figure 1. Parameters describing 
charged particles. 



the trajectories of incoming 



for r — » oo. Thus, S = +1 (S = —1) when cos a > (cos a < 
0). When cos a = 0, it is necessary to consider the second 
order term in the expansion for 9: sin 9r9 = v cos 9o (b/ro) 2 . 
Thus, in this case, S = cos 9o/\ cos 9q\. If cos$o and cos a 
both equal zero, then = 0. Note that S typically changes 
sign at points 9 where = 0. 

Since r 2 is a single- valued function of r (eq. I27p , charge 
q only reaches the grain surface if r does not reverse sign 
when r > a. From equation (|27p . r = when 



U± 



\lu 2 + Vcos9 Q + (t)' 



(29) 



where U = qQ/(mv 2 a) and V = 2qp/(mv 2 a 2 ). If the larger 
root in equation (I29p . r+, exceeds a and r > at r = r+, 
then charge q does not strike the grain. If U > 1, then charge 
q only strikes the grain if the argument of the square root in 
equation (|29p is negative, since r never reaches zero in this 
case. Assuming f > at r+, the critical impact parameter 
is given by 



(1-2(7- V cos O ) 1/2 
(-U 2 - Vcos0 o ) 1/2 



U < 1 
U > 1 



(30) 



Only trajectories with 6 < 6 cr it strike the grain surface. If 
the relevant root in equation (|30p is not real, then fe cr i t = 0. 
Note that fo cr it does not depend on the angle a. The collision 
cross section is 7r6 2 rit . 

Now we justify the assumption that r > at r 
when b = & cr i t . Differentiating equation (|18p yields 



qQ 



0. 



(31) 



The term in parentheses in equation (|3ip must vanish for all 
r, except where f = 0. Continuity implies that it vanishes at 
these locations as well, including at r = r+. Equations 
and ([31]) yield 
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if 



V cos So + (^j + U\/u 2 + V cos S + ' ''' 



at r = r + . If either U > or V cos So > 0, then clearly 
r > at r = r+, regardless of 6. If both of these quantities 
are negative, then setting b = b cr it in equation (|32[) yields 
r" = |?7 — 1| > for r = r+. 

Assuming no gas-grain drift, the mean collision cross 
section (averaged over angle So) is 

(vm 2 (1-217) , 2U + \V\<1 

a= I , 2U- \V\ > 1 (33) 

1 i7ra 2 |\/|- 1 (l-2[/ + |V r |) 2 .otherwise 

when U < 1 and 

_ i7m 2 |\/|(l - l^r 1 ^ 2 ) 2 , \V\~ l U* < 1 , , 

<0 ,\V\- 1 U 2 >1 (M) 

when U > 1. 

Integrating over the Maxwell speed distribution yields 
the factor i? from equation ((10} for a grain that does not 
drift relative to the gas: 



£(7,M) = l-7, g<2<0 and |r?| < -7 



R 



R = e~ 



R 



± {(H - 7 )(2 + H- 7 ) + 2 [l-e-^]} : 



gQ < and |7y| > —7 
sinh 1 77 



-, qQ > and ^ < 7/2 



(35) 

}. 
(36) 

(37) 



2M 

H 2 + (2- 7 )H + (2 - 7 + 7 2 /4) c _ 7/2 



I'/l 



7/2 



4|r,| 
ds I 1 



7 



gQ > and > 7 /2; (38) 

7 = qQ/(akBT gsss ), and 77 = qp/(a 2 k B T g a S ). Fig. [2] displays 
_R versus 7 for various values of \rj\. Note that equations ()35[) 
and (1371) recover the classic Spitzer (1941) expression for R 
for a charged sphere when r\ = 0. 

Since dS/dr = S/f, equations (J2TJ) and (J28J) yield 

Fx (8 ,8,VC 2 ,a) = F 2 (A,B,C) (39) 

where 

A=l + VC 2 cos O 
B = 2UC 



C = a/6 

Fi (S ,S,VC 2 ,a) 



(40) 
(41) 
(42) 
(43) 



I sin6''|d6» / S'(6'') 



^sin 2 6' - sin 2 S sin 2 a + VC 2 sin 2 S'(cos 6 - cos 0' 



(44) 



The integrand in equation (143[) is negative when 5* < 0, 
but in these cases 9 < So, so the integral remains positive. 



(32) 
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-2 



-4 - 




_l I L 



Figure 2. R vs. 7 for various values of |r)|, as indicated. 



If reaches zero at 0' = Si, then S 1 changes sign and Fi 
splits into two integrals, with limits So to Si and Si to S. 
Performing the integration in equation (|44l) . 



F 2 



F 2 = 



1 

A > 
1 

A < 0. 



VCV-B 2 +4A 



sm 



B 



VB 2 + 4A 



(45) 



:lll 



-AVC 2 -BC-A - BC-2A 



(2y^A-B)C 



(46) 



Given 2qp/(mv 2 b 2 ), 2qQ/(mv 2 b), a/b, So, and a, equa- 
tion (|39|) can be solved to efficiently find the arrival angle 
6. A less efficient, but more direct, approach is to integrate 
the equations of motion (|13p ~(|15 [) . We have written FOR- 
TRAN subroutines implementing both of these methods and 
found perfect agreement for numerous combinations of input 
parameters. 

To compute the distribution of arrival angles S for given 
values of 7 and r\, we examine a large number of trajectories 
with initial parameters So, u = v/vth, b/a, and a, where 
?Jth = (2fcB7gas/w.) 1 ^ 2 . We first select Ng values of So from 
to n, uniformly spaced in cos So. For each value of So, 
we select N v values of u, starting with u — 1.08765, the 
median value assuming the Maxwell speed distribution. We 
then select (N v — l)/2 values with it > 1.08765 spaced in 
equal-probability intervals, i.e., such that 



4 ( a 2 ( ^ 1 
—= I auu exp(— u ) = 



(47) 



Likewise for values with u < 1.08765. For each (So, it) pair, if 
fecrit > 0, then we next select Nb values of b/a between and 
fecrit /a, uniformly spaced in b 2 . Finally, for each (So,u, b/a), 
we select N a values of a, uniformly spaced between and 
2n. For each trajectory, we compute the arrival angle S. The 



6 M. E. Jordan & J. C. Weingartner 



0.8 



0.6 




0.4 - 



0.2 




Figure 3. g(cos 9) vs. cos 9 for 7 = and various values of 77, as 
indicated. 



Figure 4. g(cosS) vs. cos6> for r\ = 10 2 and various values of 7, 
as indicated. 



results are binned, with trajectories weighted in proportion 
to b 2 rit . 

Fig. [3] displays <?(cos 6), the fraction of arriving particles 
that strike with cosine of the polar angle (relative to the 
dipole moment) < cos 8, for 7 = and several values of 77. To 
construct this figure, we adopted 40 bins in 6 and No = N v = 
Nb = N a = 31. The distributions look very similar to those 
in fig. |3] when [7] < 1. Distributions for (7, —77) are identical 
to those for (7, 77), except that they are referenced to cos 9 = 
1 rather than -1. That is, (7(7, —77; cos 9) = (7(7, 77; — cos 9), 
with g(cos 9) the fraction of particles that strike with cosine 
of polar angle > cosf? when 77 < 0. As | —y | increases, the 
distribution in cos 9 becomes more uniform, as seen in Fig. 
13] for the case that 77 = 10 2 . 

Electrons arriving at the grain surface can penetrate to 
within the bulk of the grain, with an e- folding length l e ~ 
10 A (see paragraph following eq. 13 in WD01). We neglect 
this penetration since l e <C a; i.e., all arriving electrons are 
assumed to be located at r = a. 



3.3 Photoelectric Emission 

We adopt a simplified version of the procedure in WD01 
for calculating the rate at which photoelectrons are ejected 
from the grain, J pc . WD01 express the photoelectric yield 
(i.e., the probability that an electron is ejected following the 
absorption of a photon) as a product of three factors: the 
bulk yield yo, a size-dependent yield enhancement factor 771, 
and a term 1/2 that accounts for the attraction of 'attempt- 
ing' photoelectrons back to the grain when Z > 0. (Recall 
that the grain charge Q = Ze, with e the proton charge.) For 
the relatively large grains under consideration here, 771 = 1. 
The term 772 is given by (WD01, eq. 11) 



^high(^high 
1 



- 3i?k>w) / (-Eiiigh — Ei 



Z > 
z < 



,(48) 



with £i ow = —{Z + l)e 2 /a and _E high = hv - hu pct (hu pet 
is the threshold photon energy for photoemission). For sim- 
plicity, we take £"hi g h = 3 eV [assuming hv ss 11 eV and 
hfpet ~ 8 eV (WD01)], independent of Z and hv pe t- 

To find Jpc, it is necessary to integrate the photon ab- 
sorption rate _R 7 times the yield over the range of available 
photon energies above hv pc t. Since we approximate 772 to 
be independent of hv, J po oc 772. We simply choose the pro- 
portionality constant so as to reproduce the average grain 
potential of w 0.3 V from WD01. 

As for arriving electrons, we assume that holes produced 
in photemission events are located at r = a. We also neglect 
the production of an electron-hole pair when a photon is 
absorbed but a photoelectron does not escape the grain. In 
some cases, the photon absorption occurs too deep within 
the grain for the photoelectron to reach the surface or the 
photoelectron's velocity is directed away from the surface 
(resulting in 7/0 < 1)- Since l e <C a, the resulting separation 
of charge does not contribute significantly to the dipole mo- 
ment. In other photoelectron breaches the grain sur- 
face, but returns to the grain due to an attractive Coulomb 
force if Z > (resulting in 772 < 1). Such events could lead 
to a more significant change in p, but are rare; 772 ~ 0.98 
when the grain potential is 0.3 V. 



4 SIMULATIONS 

For each of the four charge transport models described in 
£|3 . 1 1 we run stochastic simulations that keep track of the 
grain dipole moment p and the orientation in space of the 
grain's rotational axis, assumed fixed with respect to the 
grain body (as would be appropriate for suprathermally ro- 
tating grains). We adopt vj_ = 1 km s _1 (resulting from 
acceleration associated with magnetohydrodynamic turbu- 



lence; Yan et al. 2004), T gas = 100 K, 



4.5 x 10" 
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s e = 0.5, T d = 15 K (hence, xo = 5 x 10~ 3 ), and y Q = 
6 x 10~ 2 (eq. 17 in WD01 with hv = 10 eV). We consider 
grains with a = 0.1 and 0.2 /im, for which the photon absorp- 
tion rate R~ t = 2.9 x 10 -2 and 5.7 x 10 -2 s -1 , respectively, 
in order to maintain the average potential at 0.3 V. We em- 
ploy a constant time step size dt, usually 31.56 s, which is 
smaller than the typical time between charging events. 

Next, we describe the simulation algorithm for perfectly 
insulating grains. At the start of each time step, we find the 
factor R (eq. I10[) for the electron collisional charging rate 
by bilinear interpolation (Press et al. 1992, p. 117) in ln7 
and In 7/, with 21 values of r\ ranging from 10 -2 to 10 2 and 
11 values of |7| ranging from 10~ 2 to 3 x 10 2 (for 23 total 
values of 7 in the tables, since both signs, as well as 7 = 
0, are included). If 7] < 10~ 2 , then a linear interpolation 
is performed in 77 alone (with 7 = 0). If 77 < 10 -2 , then 
we assume the classic Spitzer (1941) expression for R for a 
charged sphere (corresponding to r\ = 0). 

An electron arrives with probability R e dt (R e is the 
electron arrival rate; eq. I10[) . In each time-step, we check 
that R e dt < 1 (and likewise for the probability that a pho- 
toelectron is ejected). Here, as throughout the simulations, 
we use the routine ran2 from Press et al. (1992) for choosing 
random numbers. If r\ < 10 -2 , then the electron is placed at 
a random location (9, (j>) on the grain surface; 9 and <j> are the 
polar and azimuthal angles, respectively, with the rotation 
axis w = z as the polar axis. Otherwise, we interpolate to 
find the distribution function g(9'), as described in the pre- 
ceding paragraph for R; 9' is the polar angle with the dipole 
moment p as the polar axis. We choose 9' randomly from 
the distribution g(9') and the azimuthal angle <j>' is chosen 
randomly from a uniform distribution. The arrival position 
with respect to the grain body is given by 

x/a — sin #'(cos 4>' cos 9 P cos (f> p — sin </>' sin cf> p ) 

+ cos 9' sin 9 P cos <j) p (49) 

y/a = sin #'(cos 4> cos 9 V sin (j> p + sin <j> cos (f> p ) 

+ cos 9' sin 9 P sin 4>p (50) 
z/a = cos 9' cos 9 P — sin 9' cos 4>' sin d p (51) 

where 9 P and <j> p are the polar and azimuthal angles, respec- 
tively, of p relative to the grain body. 

A photoelectron is ejected with probability 7? 7 yo2/2. 
The resulting hole is located randomly on the grain surface. 

At the start of a simulation, we set p = and choose 
Z to correspond to the average potential of 0.3 V. Draine 
& Lazarian (1998) noted that a grain may have an intrinsic 
electric dipole moment due to the random orientations of 
polar constituents. Thus, our choice p = requires justifi- 
cation. 

Consider a conducting grain with N^t deep traps (model 
3 in i]3.1[) . The largest possible magnitude of the electric 
dipole moment associated with excess charges (electrons and 
holes) occupying the deep traps is p max ~ N dt ea. If the 
magnitude of the intrinsic electric dipole moment pi nt > 
Pmax, then the intrinsic dipole could not be neutralized; thus, 
flipping of p would not be possible. This situation is similar 
to that of a purely conducting grain (model 2 in jj3.1|) in 
the cold neutral medium, for which p oc Z and Z is always 
positive. 

To estimate the likely magnitude of p/ea associated 



with the intrinsic dipole, suppose each polar constituent has 

1/3 

volume Vo and dipole moment po = QeV • Assuming each 
constituent is randomly oriented, the total instrinsic mo- 
ment pint ~ N 1 ^ 2 po, with the grain volume (4/3)7ra 3 = NVo. 
Eliminating N, we find 




Even adopting relatively large values of £ ~ 0.1 and 
aV~ 1/3 ~ 500, we find pint/ea ~ 4.6, comparable (in order 
of magnitude) to the values found in the following section, 
where the instrinsic electric dipole moment pi n t is neglected. 

For any realistic grain, Ndt 3> Pint/ea; thus, we do not 
expect the intrinsic electric dipole moment to play any role 
in the long-term evolution of p, including the flipping of 
the dipole moment. Essentially, the total number of charges 
in the grain (the number of electrons plus the number of 
holes) can vastly exceed the net number of charges (number 
of electrons minus number of holes), and a slight asymmetry 
in the distribution of these charges can counter the intrinsic 
dipole moment. 

Of all the simulations with deep traps considered in 
this paper, the smallest value of Nat is ~ 100, when a — 
0.1pm and the volume per deep trap is Vt — 4 x 10 7 A . 
(This value of iVdt is almost certainly much too small to be 
realistic, but was chosen to make the computations feasible 
and to, conservatively, generate a grain with low insulating 
capability). We ran this model 3 simulation with pi nt = 
and with p lnt as estimated above; the resulting flipping times 
are nearly identical, as expected. 

At any time, the net charge and dipole moment are 
given by 

Z = N h - N e (53) 
p — Xi — Xi (54) 

i=l i=l 

where iV^ and iV e are the total number of holes and elec- 
trons, respectively, and Xi is the position of an electron or 
hole (with the origin at the grain's center of mass, i.e., the 
center of the spherical grain). In each time-step, 9 P , <j> p and 
p = \p\ are updated, if an electron arrives at or departs the 
grain. We also keep track of a ug n and a ii g n, the polar and 
azimuthal angles of the grain rotation axis with respect to 
the magnetic field direction, employing eqs. 14 and 15 from 
W06: 

d^align = fio[l — T COt 6» a li gn COs(0 a lig n + 0gyro)] dt (55) 
dflalign = — OqT Sin(0align + 0gyro) dt (56) 

where Qo is the precession rate for the case that p = (eq. 
0, 



</>gyro(t) = / dt'WgyroC*'), (57) 

Jo 

and the time-scale for gyrorotation is given by 
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Note that u; gyro varies with time, since the grain potential 
U is not constant. 

We take 0align = 0.1 and Malign = initially. Within the 
same charging simulation, we consider several different val- 
ues of uj/ujt (and thus, several different values of T; recall 
eqs. [T] [2] and |4]) . In principle, gyrorotation can affect the 
disalignment, since u! sylo fluctuates randomly as Z does so. 
However, we have found that the disalignment time is iden- 
tical for simulations that do (do not) include gyrorotation. 
Thus, we omit gyrorotation in our simulations. 

The simulations for purely conducting grains are iden- 
tical to those for purely insulating grains, except that it is 
not necessary to keep track of the electron arrival locations, 
since charge is immediately delocalized. Instead, we simply 
take p = p z z oc Z. 

For the models containing deep traps, we first specify 
the average grain volume per deep trap of a given type (i.e., 
a trap that accomodates electrons versus one that accomo- 
dates holes), Vt, then randomly place int(47ra 3 /3Vt) deep 
traps of each type throughout the grain volume. For model 
3 in H3. II (conducting grain with deep traps), an arriving elec- 
tron is immediately moved to the accomodating trap nearest 
its arrival site. This nearest trap could be a vacant electron 
trap or an occupied hole trap; in the latter case, the charges 
recombine. Likewise, the hole produced in a photoemission 
event is immediately moved to the nearest vacant hole trap 
or occupied electron trap. 

For model 4 in H3.ll (partially conducting grain with 
deep traps), each electron or hole undergoes a random walk 
through the grain, starting at its arrival location. In each 
step, the charge moves distance d rw (taken to be 30 A) 
in time i rw . Thus, for these simulations, the time step size 
dt — i r w If a charge finds itself within distance d rw of an ac- 
comodating deep trap, then it enters the trap and remains 
there until recombining when a charge with opposite sign 
arrives at the trap. The time i rw is selected as follows: 



fTcd\ 



V, 



2/3 



(59) 



with t c the typical time between charging events. With this 
choice, the typical time for a charge to travel from one trap 
to another is ~ fr a . With / ~ 1, this model lies between 
the extremes of a perfect insulator and a perfect conductor 
with deep traps. 



5 RESULTS 

We ran simulations for 11 different sets of input, with 2 re- 
alizations apiece (i.e., 2 different values of the random num- 
ber seed), for a total of 22 simulations. Table [T] displays 
input parameters for each run, as well as selected output 
parameters. We performed runs with models 1 through 3 of 
H3.ll for grain radii a — 0.1 and 0.2 /im. For model 4, only 
a — 0.1 pm is included, since the CPU time becomes pro- 
hibitive for a = 0.2 fj,m when charges execute random walks 
through the grain volume. For models 1 through 3, the total 
duration of the simulation is ttot = 10 5 yr, but substantially 
shorter t to t were obtained for model 4. 

Fig. [5] displays the component of the electric dipole mo- 
ment lying along the spin axis, p z , from a simulation for 




0.04 0.06 
t (yr) 



Figure 5. The component of the grain dipole moment (normal- 
ized to ea, the proton charge times the grain radius) lying along 
the spin axis vs. time, from a simulation of a purely insulating 
grain with a = 0.1 /im. 



a purely insulating grain with a — 0.1 (jm. Clearly, p z re- 
verses sign on a short time-scale of a fraction of a year. We 
estimate that the dipole flipping time-scale rm P ~ itot/iVfHp, 
where iV-aip is the total number of dipole flips that occur in 
the simulation. We take a flip to occur each time \p z \ in- 
creases beyond unity with p z having the opposite sign as it 
did the previous time \p z \ increased past unity. As seen in 
Table [1] Tap < 10 -3 yr for all simulations. 



The precession time flo 



2.7 x 10 -3 yr (1.1 x 



10 _ yr) for grains with a = 0.1 pm (0.2 /im). Thus, it 
is always the case that rmp < |Slo| _1 (though not always 
that Tfiip <C |f2o| _1 ). This suggests that equation 0, i.e. 

may be a good approximation for the 



T di: 



'flip 



disalignment time Tdis- Since T oc (eqs. and Q, Tdis oc 
(lo/lut) 2 - (When T < 1, as it is for suprathermally rotat- 
ing silicate grains, Tdis oc T -2 oc ui 2 when Tfn p 3> |fio| _1 as 
well; eq.|5]) As a result, extremely long simulation times are 
needed to estimate Tdis for suprathermally rotating grains. 
However, this proportionality also suggests a strategy for 
obtaining useful information with shorter simulations. We 
run for several values of ui/wt, as low as 0.1, and check 
the scaling of Tdis versus uj/lot- In fact, we find that usually 
Tdis oc (lu/lut) 2 ; thus, Tdis for highly suprathermally rotating 
grains can be estimated by extrapolation. Note, also, that 
T~ 2 |f2o| -2 is independent of the magnetic dipole moment p, 
and hence is independent of xo- Consequently, Tdi s does not 
depend on this highly uncertain parameter in this regime. 

To estimate Tdis from a simulation, we keep track of 
iVdcv, the number of times that (9 a ii gn suffers a deviation of 
1 rad. Once # a ii g n reaches a value # a ii g n,i differing by 1 rad 
from its initial value, Ndcv = L When it reaches a value 
differing by 1 rad from #align,i, Nd cv = 2, etc. When TVdov 3> 
1, Tdis ~ ttot/iVdev For the relatively small values of lu/ujt 
under consideration, this condition obtains. However, this 
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Table 1. Simulation Parameters and Outputs 



Modefi qp Vt^_ Run \p z Uj_ r fliB £, hot*. l.q* 1.5 4 2.0* 

(im A 3 s ea 10 -4 yr yr 



1 


0.1 






1 


2.21 


5.2 


1.0E5 


2348 


218 


30 


1 


0.1 






2 


2.18 


5.2 


1.0E5 


2418 


195 


26 


2 


0.1 






1 


2.11 




1.0E5 


72 


5 


1 


2 


0.1 






2 


2.12 




1.0E5 


44 


6 


2 


3 


0.1 


4.0E7 




1 


1.77 


6.8 


1.0E5 


1822 


179 


19 


3 


0.1 


4.0E7 




2 


1.77 


6.2 


1.0E5 


1800 


160 


19 


3 


0.1 


4.0E6 




1 


2.11 


5.8 


1.0E5 


2451 


220 


26 


3 


0.1 


4.0E6 




2 


2.13 


5.8 


1.0E5 


2358 


228 


32 


1 


0.1 


4.0E7 


2.3 


1 


2.16 


2.2 


7.29E2 


14 


1 





4 


0.1 


4.0E7 


2.3 


2 


2.35 


2.5 


7.29E2 


22 


1 





4 


0.1 


4.0E7 


1.15 


1 


2.17 


2.5 


3.64E2 


6 


1 





1 


0.1 


4.0E7 


1.15 


2 


2.09 


2.4 


3.64E2 


11 


1 





1 


0.1 


4.0E7 


9.2 


1 


2.56 


2.4 


2.92E3 


81 


11 


1 


4 


0.1 


4.0E7 


9.2 


2 


2.56 


2.5 


2.92E3 


76 


7 


1 


4 


0.1 


4.0E6 


2.3 


1 


2.07 


5.0 


7.29E2 


17 


2 





1 


0.1 


4.0E6 


2.3 


2 


2.11 


5.1 


7.29E2 


19 


3 





1 


0.2 






1 


2.23 


2.6 


1.0E5 


172 


16 


2 


1 


0.2 






2 


2.18 


2.6 


1.0E5 


150 


12 


1 


2 


0.2 






1 


2.00 




1.0E5 


4 








2 


0.2 






2 


2.00 




1.0E5 


4 


1 





3 


0.2 


4.0E7 




1 


2.07 


2.9 


1.0E5 


164 


17 


2 


3 


0.2 


4.0E7 




2 


2.07 


2.9 


1.0E5 


139 


14 


2 



a Number of 1 rad deviations in alignment angle a ii gn . 

^ Prom section [XT] 

c Grain radius. 

^ Volume per deep trap. 

e Duration of random walk step. 

f Average of the absolute value of the component of the electric dipole moment lying along the spin axis 
(normalized to ea, the product of the proton charge and the grain radius). 
9 Estimate of the electric dipole moment flipping time. 
^ Duration of the simulation. 
1 Suprathermality u>/u)j-. 



is not typically the case when uj/ujt^W- In general, we 
estimate r d i s ~ t to t/(Nd cv + |A# a i ign |), where A# a i ign is the 
value of f?aii gn at the end of the simulation minus its value at 
the last time iVdcv was incremented (which may have been 
the start of the simulation, if iVdev = 0). Table [T] indicates 
the values of TVdev for log 10 (uj/LJT) = 1-0, 1.5, and 2.0. Of 
course, the resulting estimate of Tdi s is not very reliable for 
the cases where iVdcv ~ 1. 

For each simulation, we keep track of # a ii gn for 7 dif- 
ferent values of log 10 (u;/t<;T), evenly spaced between -1.0 
and 2.0. Fig. [6] displays cos^aijgn versus time from a sim- 
ulation of a perfectly insulating grain with a — 0.1 /im and 
ui/lut = 10 2 . 

Fig. [7] displays Tdi s versus lu/ujt for perfectly insulat- 
ing grains with a = 0.1 and 0.2 ^m. For each case, Tdis is 
taken to be its average over the 2 realizations. The solid 
(dashed) curves are Tdi s from equation (0 for a = 0.1 ^m 
(0.2 fim). We employ raip and the average value of \p z \ (for 
use in evaluating T) as determined from the simulation. The 
agreement between the measured values of Tdi s and those 
calculated with equation is surprisingly good. The ex- 
pectation that rdi s oc (uj/lot) 2 is well confirmed. 

The disalignment times found using equation and 



the values of rflip from the simulations are substantially 
shorter than those from W06 (see figs. 2 and 3 in W06); the 
discrepency exceeds 2 orders of magnitude when a — 0.1 /jm. 
Our simulations yield much larger values of |p z | av than es- 
timated by W06, and rdis oc Ipzlav 2 . The estimate of rm P 
in W06 is also substantially larger than our result. When 
equation is used (rdi s oc T~j p ), this partially compensates 
for the difference associated with the p z estimates. However, 
given the larger estimate for Ts.i p , W06 employed equation 
when a — 0.1 /jm; in this case, Tdi s is larger by a factor 
~ 2 when equation ((8} is used than when equation is 
used. 

Fig. [S] shows the ratio of Tdi s for several simulation runs 
to its value for the perfectly insulating case, t 111s , for a — 
0.1 fim. All of the simulations from Table Q] are included, 
except for the perfectly conducting case. The results always 
lie within w 50 percent of unity, with somewhat greater 
scatter when uj/ujt > 10; the results for these high-w cases 
are not particularly reliable, since the corresponding iVd ev 
are small (see Table 0. The ratio rdis /Tins also lies within 
50 percent of unity for the model 3 run for a = 0.2 fj,m. 

Fig. [9] displays the ratio of Tdis for a perfectly conduct- 
ing grain, T con d, to Ti ns , for a = 0.1 and 0.2 /im. We assumed 
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0.2 0.4 0.6 0.8 1 
t (10 5 yr) 

Figure 6. cos 9 a i ign versus time from a simulation of a perfectly 
insulating grain with a = 0.1 /im and uj/ujt = 10 2 ; N^ cv = 30 for 
this case. 
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Figure 8. Ratio of disalignment time T<j; s for various models to 
that for a perfectly insulating grain, Ti ns , for a = 0.1 ^m. Model 
number from SI3.1I volume per deep trap, Vt\ and random walk 
time, t rw , are indicated. 
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Figure 7. Disalignment time vs. suprathermality for purely in- 
sulating grains. 



60 



>40 - 



20 



_L 



, — i — | — i — i— i— i— |- 
a=0.2/um 

A a = 0.1/i.m 



-10 12 
l°gio("/ w T) 

Figure 9. Ratio of disalignment time for a perfectly conducting 
grain, r con( j, to that for a perfectly insulating grain, T; ns , for a = 
0.1 and 0.2 /^m, as indicated. 



that p z =0.1 Zea, which seems conservative for grains suf- 
ficiently asymmetric to produce the observed polarization. 
However, a solution of the electrostatic boundary value prob- 
lem for model aspherical grains would be needed to confirm 
this choice. The disalignment times tend to be 1 to 2 orders 
of magnitude longer for conducting grains than for insulating 
grains. This is not surprising, since Z, and hence p z for con- 
ductors, does not change sign (although it does fluctuate). 



Note that the data points for uj/uit^W are not reliable, 
given the small values of JVdev in these cases (Table H}. 



6 CONCLUSIONS 

We have conducted a more detailed analysis of grain dis- 
alignment associated with the time-varying electric dipole 
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moment than was attempted in W06, focusing on suprather- 
mally rotating silicate grains. We considered 4 idealized 
models for how charge is transported within the grain f £!3,l[) : 
a perfect insulator, 2 models involving special sites in the 
grain ('deep traps') where electrons or holes are effectively 
trapped, and a perfect conductor. The resulting disalign- 
ment times Tdis for the first 3 models are highly consistent 
(Fig. [8} and substantially shorter (up to 2 orders of magni- 
tude) than those obtained by W06 (cf. Fig.[7]here with figs. 
2 and 3 in W06). We expect the behavior of real grains to 
be bracketed by these 3 models. Disalignment proceeds more 
slowly (up to 2 orders of magnitude; Fig. [9]) for conducting 
grains, but we do not expect this idealization to be realistic 
for interstellar grains. 

In treating the collisional charging, we neglected the 
gas-grain drift. Drift can, in principle, affect the time varia- 
tion of the electric dipole moment. For a non-rotating grain, 
there may be a stable contribution to p directed along the 
drift velocity. For a grain rotating uniformly about cii, the 
charging rate may have some dependence on latitude on the 
grain, suppressing flips in p z . We examine this possibility in 
Appendix lAl and conclude that flipping is not suppressed. 

In the radiative torque alignment scenario, suprather- 
mal rotation with to /lot ~ 100 and alignment times > 10 5 yr 
appear to be typical (Draine & Weingartner 1997; Lazarian 
& Hoang 2007; Hoang & Lazarian 2008), though additional 
studies are needed to confirm these results. We have found 
disalignment times < 10 J yr when lo I lot ~ 100 (Figs. [7] 
and [8]), presenting a severe challenge to the radiative torque 
model. 

Much of the physics involved in the disalignment mech- 
anism has not been directly verified, including the details of 
the charging and the turbulence-induced grain acceleration 
(Yan et al. 2004). Perhaps current models of these processes 
are incomplete in such a way as to overestimate the magni- 
tude of the disalignment. 

Alternatively, interstellar grains might contain super- 
paramagnetic inclusions (Jones & Spitzer 1967) , which could 
increase the magnetic susceptibility by orders of magnitude. 
The parameter T would be decreased by the same factor, 
and the disalignment time oc Y~ 2 when T 1 (eqs. [7] and 
[8|. Recently, Lazarian & Hoang (2008) found that the pres- 
ence of superparamagnetic inclusions can modify alignment 
by radiative torques, yielding a higher degree of alignment 
than experienced by grains free of inclusions. Perhaps su- 
perparamagnetic inclusions also suppress drift-induced dis- 
alignment. 
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APPENDIX A: COLLISIONAL CHARGING 
FOR A DRIFTING GRAIN 

Consider a grain drifting with velocity v gr with respect to 
the gas. The grain rotates uniformly about &i, which is in- 
clined at angle 9j v relative to v sr . 



To treat the collisional charging in this case, we first 
construct a large sphere with radius rt,i g instantaneously 
centered on the grain. Adopting the rest frame of the gas 
and taking the direction of the drift velocity u gr as the polar 
axis for spherical coordinates, the velocity v of a gas-phase 
particle has components (v, 9[ n , (j>in)- The rate at which gas- 
phase particles enter the large sphere from within solid angle 
dcos#i n d0in about (#i n ,</>m) and with speeds between v and 
v + dv is 



dR = ivr bi „n — d cos 9 ln d(f>i n P(v)dv \v — v g 



(Al) 



where P(v) is the Maxwell speed distribution. After inte- 
grating over (/>i n , 



2 /8fc B T\ 1/2 
an = nr bis n I 



V 7rm / 



dR 



with 



dR = du d cos #i n u 2 u\ exp(— u 2 ); 

the dimensionless speed u = v/vth and 

Itl = (v 2 + U 2 + 2-UMgr COS Sin) 1 



(A2) 



(A3) 



(A4) 



is the particle's dimensionless speed in the rest frame of 
the grain (u gI — v SI /v t h)- Integrating over the entire large 
sphere yields 



R= 1 + 



u k [ diexp(-u gr x)(l-^) 3 . (A5) 



The arrival angle 9\ in the rest frame of the grain is charac- 
terized by 



costfi = 



u cos #i n + it g 



(A6) 



When v± = 1 km s" 1 and T gas = 100 K, u gr = 0.0182 
for electrons and 0.778 for protons. For each of these values 
of M gr , we calculate ui and dR for 10 10 values of (u, 9i n ) (10 s 
for each input variable) spaced evenly in probability (as de- 
scribed in the text surrounding eg . I47fl and with a maximum 
value of u = 4. The results are used to construct the prob- 
ability P(u\) that an incoming particle has dimensionless 
speed in the grain's frame < u\ (with 100 bins in ui). For 
each value of u\, the corresponding cumulative probability 
P(cos#i) is constructed, again with 100 bins. Note that the 
minimum possible value of cos#i is —1 when u > u SI and 
[1 — (u/ugr) 2 ] 1 ^ 2 when u < it gr . 

To simulate the collisional charging, we adopt a time 
step 10 times smaller than the inverse of the rate at 
which electrons enter the large sphere surrounding the grain 
(eqs. IA2I and IA5[1 . In each step, we draw a random num- 
ber to determine whether or not an electron enters the large 
sphere; likewise for a proton. When a charged particle en- 
ters the large sphere, a value of ui is picked randomly from 
its distribution. Then, cos 9i is chosen randomly from the 
distribution for the given u\. The final component of the 
particle's velocity in the grain's rest frame, cj>i, is selected 
randomly from a uniform distribution between and 2tv. 
Two components of the incoming particle's position remain 
to be determined (given r = rbi g ): the impact parameter b 
and the azimuthal angle ai . These are both chosen randomly 
(b from a uniform distribution in ft 2 ). 

At this point, the position and velocity of the incoming 
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particle are specified relative to a coordinate system at rest 
with respect to the grain and with v SI as the polar axis ('i> g r- 
coordinates') . Denoting Cartesian axes in this coordinate 
system as x v , y„ , and z„ , the velocity is given by 

«i = — vtWU\(x v sin 0i cos0i+y„ sin 9\ sin 4>i+z v cos#i)(A7) 

and the position by 

x v = 6 cos Qi cos 9i cos </>i — 6 sin qi sin </>i 

+Zarr sin 6>i COS 0i (A8) 

y v — b cos ai cos 81 sin 0i + b sin ai cos 0i 

+z arr sin #1 sin 0i (A9) 

z v — ~b cos ai sin 9i + z arr cos #1 ( A10) 

where z arr = (rl ig - fe 2 ) 1/2 . 

Denoting a Cartesian coordinate system attached to the 
grain body by (xj,y.j,zj), 

zj = 01 = z v cos 9 Jv + x v sin 9 Jv (All) 

4j = (£„ cos 6j v — Zv sin 9j v ) cos A + y„ sin A (A12) 

yj = j/u cos A — (x v cos #j„ — z„ sin #j„) sin A (A13) 

where A is the phase angle of the grain's rotation and is se- 
lected randomly. (We neglect the rotation of the grain during 
the approach of the gas-phase particle.) Cartesian axes with 
p as the polar axis ('p-coordinates') are given by 

z p = p = xj sin 9 pJ cos </> p j+yj sin 9 pJ sin 4> p j+zj cos pJ (A14) 

xp = xj cos9 pJ cos4> pJ +yj cos9 p j sm4>pj—zj sm8 p j(Al5) 

yp = -4 j sin 4>p.j + yj cos (A16) 



dr 1 / dx p dj/p dz p 



From equations (IA1 A16|) . we find the following dot prod- 
ucts for use in transforming the position and velocity of the 
incoming gas-phase particle from i> gr -coordinates (eqs. IA7l - 
IA10|) to p-coordinates: 



Xp 


' Xd 


= cos 9 p j cos 9j v cos(<t> p j + A) — 


sin 


9 p j sin 


9 Jv (A17) 


Xp 


• yv 


= cos 9 p j sin(0 p j + A) 






(A18) 


Xp 


■ Zv 


= —cos9 p j sm9.j v cos(4>pj+X)- 


sin 


9 p j cos 


e Jv (Ai9) 


Vp 


Xv 


= — cos 9j v sm((j)pj + A) 






(A20) 


Vp 


yv 


= cos((p p j + A) 






(A21) 


Vp 


Zv 


= sin 9j v sm(4> p j + A) 






(A22) 


Zp 


Xv 


= sin 8pj cos 9j v cos(0 p j + A) + 


cos 


9 P j sin ( 


hv (A23) 


Zp 


Vv 


= sin 9pj sin(0 p j + A) 






(A24) 


Zp 


Zfj 


= — sin 9 pJ sin 9 Jv cos(0 p j+A) + 


cos t 


) pJ cos t 


hv{A2h) 



Finally, the Cartesian p-coordinates of the incoming 
particle's position and velocity are used in the following ge- 
ometric relations to find the components in spherical coor- 
dinates: 



9 — cos _1 (z p /r) 



2 tan 



TpO Xp 

Vp 



(A26) 
(A27) 

(A28) 



d9 
dt 

dcf> 
lit 



1 

rr p o 



(x P 



dzp 
dt 

dyp 
dt 



dr\ 

- Zp Tt) 



dxp 
dt 



(A29) 
(A30) 
(A31) 



with r p o = (Xp + y^) 1 . The critical impact parameter 6 cr i t 
depends on the particle speed Vi = WttWi and cos^o = 
— v^ 1 dz p /dt (eq. l30p . lib < 6 C rit, then we integrate the equa- 
tions of motion (|13[1 - (|15[) to determine where on the grain 
surface the particle hits. 

We have tried various values of rbi g . Of course, larger 
values yield higher accuracy but also require smaller time 
steps. We found that r^g = 50a yields high accuracy and is 
not prohibitively time consuming. 

Substituting the collisional charging procedure de- 
scribed here in our charging simulations (and including both 
electrons and protons), we examined a perfectly insulating 
grain with a = 0.1 /im. With a duration of 100 yr, we found 
that |p z |av/ea w 2.5 and Tflj p ranges from 5.8 x 10 -4 to 
6.0 x 10~ 4 as cos9j v ranges from to 1. These are very 
close to the results obtained previously, neglecting grain drift 
(and ignoring protons) in the treatment of collisional charg- 
ing (see Table [TJ. Due to precession of J about B, the angle 
9j v changes on a time-scale short compared with the sim- 
ulation time of 100 yr (but an order of magnitude longer 
than Tflip). However, our results imply that the behavior of 
the electric dipole moment is insensitive to the value of 9j v . 
Thus, we conclude that the neglect of grain drift in ^5] does 
not yield significant error. 
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